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We examine the optical conductivity in antiferromagnetic (AFM) iron pnictides by mean- 
field calculation in a five-band Hubbard model. The calculated spectra are well consistent with 
the in-plane anisotropy observed in the measurements, where the optical conductivity along the 
direction with the AFM alignment of neighboring spins is larger than that along the ferromag- 
netic (FM) direction in the low-energy region; however, that along the FM direction becomes 
larger in the higher-energy region. The difference between the two directions is explained by 
taking account of orbital characters in both occupied and unoccupied states as well as of the 
nature of Dirac-type linear dispersions near the Fermi level. 
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Iron pnictide superconductors"'^^ have recently at- 
tracted much attention and have been providing a new 
direction in the studies of high-temperature supercon- 
ductivity. One significant property of iron pnictides is 
the existence of the antiferromagnetic (AFM) metaUic 
phase in their parent compounds. This imphes a close 
connection between superconductivity and magnetism in 
high-temperature superconductors. In order to elucidate 
the mechanism of high-temperature superconductivity in 
iron pnictides, an understanding of the properties of the 
AFM metallic phase is strongly desired. 

Recent resistivity measurements of detwinned sam- 
ples of the AFM phase^""*) in the 122 systems, AFe2As2 
(A=Ca, Sr, Ba), have demonstrated anisotropic electric 
transport in the Fe-As layer: Below an orthorhombic 
structural transition temperature, the electric conduc- 
tivity along the direction with the AFM alignment of 
neighboring spins (the x-direction in our notation below) 
is larger than that along the ferromagnetic (FM) direc- 
tion (y-direction). Among the above three compounds, 
BaFe2As2 has revealed an anisotropy even in a tetrago- 
nal phase above the structural transition. This suggests 
the presence of a nematic order prior to the structure 
transition and magnetic ordering, which might be related 
to the multiorbital nature of iron pnictides. 

In addition to the electric resistivity, the optical con- 
ductivity of the detwinned BaFe2As2 has shown an 
anisotropic spectral behavior:^' The optical conduc- 
tivity along the x-direction is larger than that along 
the y-direction below 180 meV while that along the y- 
direction becomes larger above 180 meV. Together with 
recent scanning tunneling microscopy measurements,^-* 
the study of electronic excitations has revealed a large 
electronic anisotropy in the AFM metallic phase. There- 
fore, it is crucial to clarify the origin of the electronic 
anisotropy, focusing on the orbital nature of electronic 



excitations coming from Fe3(i bands. 

In a theoretical work on optical conductivity by the 
dynamical mean-field theory (DMFT) combined with 
first-principles band-structure calculation, Yin et alJ^^ 
have reported that the anisotropy in electric transport 
originates from Pauli's exclusion principle, which blocks 
the transition to the FM direction. However, their as- 
sertion seems to be too simple to give a full explana- 
tion of the anisotropy in the entire energy range. A 
very recent fixed-spin-moment band-structure calcula- 
tion has yielded anisotropic optical conductivity quali- 
tatively consistent with experimental data.^^ However, 
the microscopic origin of the anisotropy has not been 
fully discussed. Therefore, no clear explanation has been 
proposed yet. To solve this problem, we provide here 
an interpretation that brings a new viewpoint of the 
anisotropy based on our calculations. 

In this study, we investigate the anisotropy of the opti- 
cal conductivity by mean-field calculation in a five-band 
Hubbard model. Although we use crude approximation 
as compared with the DMFT work,^^ physical interpreta- 
tion is easily made to examine the origin of the anisotropy 
of the optical conductivity. Our calculations are consis- 
tent with the anisotropic behavior observed experimen- 
tally. Taking account of the symmetry of the orbital 
character and the presence of Dirac-type linear disper- 
sions, we obtain a physically transparent origin of the 
anisotropy, which depends on the energy region of the 
optical conductivity. We also find that AFM ordering is 
sufficient to give adequate anisotropy in the optical con- 
ductivity, and thus orbital order is not necessary at least 
below the Neel temperature. 

Considering an Fe square lattice, we start with a 
multiband Hubbard Hamiltonian for a d-electron system 
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H = Ho + Hi. Here, 



Hn = 



A 



(1) 
cre- 



is the five-band hopping Hamiltonian, where 
ates an electron with a wave vector k and a spin a at 
an orbital /i, and A — (A^^, Aj,). t(Aa;, A^; /z, i') and 
are the in-plain hopping integrals and on-site energies, 
respectively, presented by Kuroki et al.^^ Hi is the inter- 
action Hamiltonian expressed as^*^) 



Ht^U 



{U-2J) J2 



U-3J 



Et t 
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where i is the Fe-site index, U is the intraorbital Coulomb 
interaction, J is the Hund coupling, and the pair hopping 
is set to J. The Fe-Fe bond length is set to unity and the 
a;- and y-directions are along to the nearest Fe-Fe bonds. 
Our calculation uses only two dimensional hopping in- 
tegrals. Since interplane hopping integrals are estimated 
to be small, this would hardly affect the optical con- 
ductivity discussed below. 

We self-consistently solve mean-field equa- 
tions with the order parameter defined by 
Cfc^_Q^o.Cfcyo.) with the order- 
ing vector Q and the number of k points, N, in the first 
Brillouin zone (BZ) of the paramagnetic phase. For the 
observed AFM ordering, we take Q = (tt, 0), implying 
the AFM arrangement of spins along the x-direction 
and the FM arrangement along the y-direction. The 
very small difference of the bond length between the 
X- and y-directions in the orthorhombic phase is not 
taken into account for simplicity in our calculation.^^) 
In addition to Q = (tt, 0), the vector Q = (0, 0) that can 
represent orbital ordering^^^ is included. The quasipar- 
ticle state 7^^^ = I;q=(o,o),(^.o) '^P^.eaik + Q)c1^q^^ 
diagonalizes the Hamiltonian with the eigenvalue -Efceo-, 
where e is the band index. The average (• • • } is taken at 
absolute zero in our calculation. The computations are 
performed on a system with N — 500 x 500. 

We set U = 1.2 eV and J = 0.23 eV to yield a 
magnetic moment M = 0.85/iB (Bohr magneton) close 
to the measured moment of the 122 system, where 

^ = I];.("(^,o)MMt ~ "-(Tr,o)MMi)MB- Wc uote that orbital 
off-diagonal order parameters result in zero values. Or- 
bital ordering may be measured using Uyz — Uzx with 
— X]<j(^*(o.o)aia'ct): whose value is found to be small 
(-0.062). 

Interband contributions to the real part of the optical 




Fig. 1. (Color online) Interband contribution to the optical con- 
ductivity along the x-direction a^xi}^) (red solid line) and in- 
direction ayy{u)) (blue dashed line). The Drude component is 
not included. The energy scale of the transverse axis is obtained 
by multiplying the original scale in the mean-field calculation by 
1/3 in order to simulate band renormalization due to the cor- 
relation effect. The interband contributions are separated into 
four regimes labeled as A: range of 0-10 meV, B: 10-65 meV, C: 
65-147 meV, and D: 147-200 meV. 



conductivity are expressed as""^^) 
— TT / e \ 2 

fe.( 



aap (w > 0) 



E [f{Ekea)-f{Eke'.)\ 



(/3) 
ee'crJ 
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where e is the elementary charge, / is the Fermi dis- 
tribution function, and Ckll'a arising from the current 
operator has the form 

x^^e.(fc + Q)C.'.(fc + Q), (4) 

with A*^") the a component of the vector A. The Drude 
component coming from the intraband transition can be 
obtained in the limit w -> and written as 



CTa/3(w = 0) = 



1 



E 

— Alp ,e,cr 



W 1 

eecrJ 



vik) 
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where Np is the number of k points on the Fermi surface 
and v{k) = \dEkea/dk\. 

Figure 1 shows interband contributions to the opti- 
cal conductivity. In the figure, we use a renormalized 
energy scale with a factor 1/3 to the energy axis, i.e., 
0.1 eV in the figure corresponds to 0.3 eV on the original 
scale of energy. This factor is taken from comparisons 
between the dispersion observed by angle-resolved pho- 
toemission and the theoretical dispersion determined by 
first-principles calculation. The factor corresponds 
to a band renormalization effect that is included in nei- 
ther the first-principles calculation nor our mean-field 
calculation. The fact that the scaled result in Fig. 1 is 
consistent with the experimental data (shown below) in- 
dicates a crucial role of the self-energy due to the corre- 
lation effect. 

In Fig. 1, there are characteristic structures separated 
into four regimes. In the range of 0-10 meV (label A), 
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Fig. 2. (Color online) Distribution map of interband transitions 
contributing to the optical conductivity for each regime (A-D) 
in the momentum space: A (red +), B (green x), C (blue *) and 
D (purple □). The momentum positions giving rise to dominant 
contributions for each regime are denoted by A to D. 
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Fig. 3. (Color online) Band dispersions of the Q = (tt, 0) AF 
ordered phases in the five-band Hubbard model. The labels A, 
C, and D represent the direct excitation process predominantly 
contributing to the optical conductivity in regimes A, C, and D, 
respectively. The label B' denotes the position contributing ayy 
in regime B. 



ayy is small (close to zero) and shows no characteristic 
excitation peak, while cr^x shows a peak. In the range 
of 10-65 meV (label B), ayy increases more slowly than 
axx does. In addition, there is a plateau in axx- In the 
range of 65-147 meV (label C), a^x increases, and a peak 
appears at 80 meV. In the range of 147-200 meV (label 
D), ayy becomes larger than axx in contrast to that in 
the lower-energy regimes. 

Since the optical conductivity is given by momentum- 
conserved interband transitions, it is easy to assign mo- 
mentum regions contributing to regimes of A-D. Fig- 
ure 2 shows the regions where the energy of direct in- 
terband transition across the Fermi level is equal to the 
energy range of each regime. Note that the distribution 
in the momentum space is symmetric with respect to 
kx = 7''/2 because of the Q = (tt, 0) ordering. Regime A 
is concentrated in a very narrow region close to the ky 
axis. Regime B is distributed around k = (0.27r,0) and 
(O.Stt, 0). These momenta are close to a crossing point 
of Dirac-type linear dispersions. ^^"^^^ Regime C is as- 
signed to momenta around (0,0) and (0,7r), and regime 
D surrounds regime C. The points denoted by A to D 
represent the momenta making a dominant contribution 
to the largest spectral intensity within each regime, i.e., 
axx at a; = 5 meV for regime A, axx at a; = 25 meV for 
regime B, axx at = 85 meV for regime C, and ayy at 
w = 160 meV for regime D. We find that the positions 
are located along the highly symmetric line in BZ, except 
for B. This fact is a key to determining the difference be- 
tween axxi^) and ayy{uj), as discussed below. 

In order to clarify the k points contributing to the op- 
tical conductivity, we illustrate band dispersions along 
the symmetric lines in Fig. 3 together with their posi- 
tions. The label B' represents an excitation contributing 
to a peak of ayy at = 63 meV in regime B (see Fig. 1). 

On the symmetric lines in BZ, it is convenient to clas- 
sify wave functions of the initial and final states in terms 
of the symmetry of the lines. The crystal structure of 
iron pnictides has a C2 rotation about the x- and y-axes. 
When the sign of the orbital wave function changes (re- 
mains) under the C2 rotation, we define the parity of 
the function as — (+). For example, the parity of the 



yz-orbital about the y-axis is — . We analyze the wave 
functions of the initial and final states, and summarize 
their parity and orbital components in Table I. The ini- 
tial and final states are connected by a dipole operator 
associated with the applied electric field. The parities of 
the electric field about the a;-axis rotation are -I- and — 
for the X- and y-polarizations, respectively, and vice versa 
about the y-axis rotation. Combining the parities of ini- 
tial and final states with the parity of the electric field, 
we can judge whether or not the interband transition at 
A to D is possible. 

Regime A. — axx mainly comes from A on the (7r,0)- 
(7r,7r) line parallel to the y-axis in Fig. 3. The orbital of 
the initial (final) state has the parity -I- (— ) about the 
C2 rotation around the y-axis. The parities of the electric 
field along the x- and y-axes are — and +, respectively. 
Therefore, ayy at A vanishes because -f (initial) times -|- 
(dipole) is not — (final), while axx at A is allowed because 
+ (initial) times — (dipole) is equal to — (final). This is 
the main reason why axx is larger than ayy in regime 
A. In this regime, we cannot compare our result with 
experiments where the Drude component is dominant. 
The Drude component is discussed below. 

Regime B. — The optical conductivity mainly comes 
from the momentum space around B in Fig. 2. The spec- 
tral behavior in this regime is explained well by the pres- 
ence of the Dirac dispersions near B. In the vicinity of the 
Dirac point, the velocity is constant because of linear dis- 
persion, leading to an energy-independent spectral inten- 
sity. This is consistent with the appearance of a plateau 
in axx- Note that, since the Dirac point is at —10 meV 
from the Fermi level in the (0,0)-(7r,0) direction, the 
contribution of this Dirac cone emerges above 10 meV. 
In the experimental data, plateaus appear from 45 to 
75 meV in both axx and ayy.'^'' Our calculations show a 
plateau in axx , but not in ayy . This also emerges in a pre- 
vious DMFT work.^^ In our mean-field calculation, the 
velocity of the Dirac cone along the x-direction is nearly 
three times as large as that along the y-direction.^'^) Such 
a difference in the velocity contributes to the difference 
between axx and ayy. We should note that the Dirac 
cones observed in experiments^^-* seem to be not very dis- 
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Table I. Orbitals of initial and final states in each position. The sign in the brackets represents the parity of C2 operation around the 
rotation axis. 



position 


rotation axis 


initial 


final 


A 


y 


<^322 -r2 , dzx ,d^2_y2 (+) 


dyz, dxy{-~) 


B' 


X 


dzx 1 dxy ( ) 


d^z^ -r2 , dy z , d^2 _y2{-\-) 


C 


X 


dzx 1 dxy ( ) 


dzx 1 dxy ( ) 


D 


y 


dyz, dxy{—) 


dyzi dxy{~) 



torted, as calculated.^*) Reflecting this, the measured op- 
tical conductivity shows less anisotropy in this energy re- 
gion unlike that in our calculations. The interband transi- 
tions along the (0,0)-(7r,0) direction in Fig. 3 are allowed 
for (Tyy, as shown in Table I. As a result, the transition 
B' appears as a peak at approximately w = 65 meV in 
(jyy, as shown in Fig. 1. 

Regime C. — At C in Fig. 3, the wave hmctions of both 
the initial and final states have the — parity about the 
X-axis rotation, as listed in Table I. Both states can be 
connected by the x-polarized electric field, but not by 
the y-polarized electric field. This explains why a^x is 
larger than ayy. Unlike in the experiments, peak struc- 
tures emerge at 110 mcV in a^x and at 125 meV in 
cTyy^^ These peaks correspond to a sharp peak at ap- 
proximately 85 meV in Gxx and a broad peak at approx- 
imately 95 meV in (jyy in our calciilations. 

Regime D. — The main contribution to Oyy comes 
from D on the axes parallel to the y-direction in Fig. 3. 
This is similarly explained using the symmetry argument 
that the initial and final states are connected by the 
^/-polarized electric field. The opposite anisotropy com- 
pared with those of regimes A and C is, thus, caused by 
different nature of orbitals in the initial and final states. 
Our results are consistent with the experimental data 
showing that (Tyy is larger than a^x iu this regime. 

The Drude component coming from intraband transi- 
tion is calculated from the present band structure. The 
calculation leads to a larger Drude component along the 
x-direction with the ratio Oxxi'-^ = 0)/ayy(a; = 0) = 1.3. 
This difference can be understood from the anisotropy of 
the Fermi velocity.^") 

In order to extract the effect of orbital ordering on 
the anisotropy of the optical conductivity, we performed 
a calculation without the Q = (0,0) order parameter 
(not shown). We find that the anisotropy remains almost 
unchanged, and thus conclude that orbital ordering is 
irrelevant as the origin of the anisotropy but that AFM 
ordering plays an essential role as the origin. 

In summary, we have investigated the anisotropy of the 
optical conductivity by mean-field calculation in a five- 
band Hubbard model. Our results are consistent with the 
anisotropic behavior observed experimentally. The origin 
of the anisotropy is explained in terms of the orbital char- 
acters of the initial and final states contributing to inter- 
band transitions. In addition to the orbital symmetry, 
the presence of Dirac-type linear dispersions is crucial 
for understanding a plateau-like spectral shape in the 
optical conductivity. We have found that AFM ordering 
is enough to give adequate anisotropy in the optical con- 



ductivity, and thus an orbital order is not necessary at 
least below the Neel temperature. The fact that the or- 
bital degree of freedom plays a crucial rule in the AFM 
metallic phase suggests that the orbital character should 
be taken into account in further study of superconduc- 
tivity in pnictides. 
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